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ABSTRACT 


Oceanographic and meteorologic data, namely expendable 
or mechanical bathythermograph traces and wind reports have 
been systematically collected for years and every day more 
observations are made. 

To handle such a large amount of information, it is 
advantageous to set up a computerized model which is: (1) 
automatic, (2) performs an objective analysis of XBT/BT 
traces to identify the most significant thermal features, 
Such as thermoclines and inversions, (3) restricts the human 
participation to an initial stage dealing with organization 
OL Gata available and. selection of proper values for the 
parameters of the moaei, and (4) preserves the initial 
information. 

A FORTRAN IV program takes into consideration the above 
requirements to: (1) perform objective analyses of XBT/BT 
digital traces, (2) process wind information, (3) study the 
statistical characteristics of the original and processed 
data, and (4) plot the generated fields. 

The digital model was tested and its application may be 
extended to the study of any environmental variable continu- 


ously distributed along a vertical or horizontal axis. 
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Lo --OBJECTIVE..OF THE. THES LS. 


The study of oceanographic phenomena occurring in an 
area west of the Strait of Gibraltar from Cabo Finisterra 
(in Europe) to Canary Islands (in Africa) (Figure 1) based 
on Climatological information, constituted the initial pur- 
pose .of<this .thesis. 

Later, the availability of oceanographic and meteorologic 
data from Fleet Numerical Weather Central (in Monterey), 
broadened the descriptive nature of the above intention 
Siveng Opportunity to investigate the variations™o£f—the 
oceanic thermal structure taking place in the upper layers. 

The large amount of accessible data, namely expendable 
bathythermograph traces and surface uncorrected geostrophic 
wind fields, were processed using a digital model and pre- 
sented in a form simple to analyze. Initially, the model 
was set up to study the distribution of temperature with 
depth; soon, it became evident that it could be applied to 
other studies dealing with the continuous vertical distribu- 
tion of any environmental parameter (sound velocity, salinity, 
laght extinction, etc.). 

The lack of information about air-sea temperature differ- 
ences made impossible the transformation of the uncorrected 
geostrophic wind into "friction velocity" and cut short the 
last part of the study, which was to establish the statistical 
interrelation between the corrected wind information and the 


thermal features given by the objective analy<is of traces. 
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Figure 1. Eastern North Atlantic Ocean 
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The new objective of the thesis then is defined in two 
parts: 

a. To develop a digital model to analyze the continuous 
vertical distribution of any environmental parameter, to 
process the wind data and to present the generated informa- 
tion in a tractable form to reveal geophysical phenomena 
and their statistical connections. 

ba To test the capability of the computer program by 


means of data available. 
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II. OBJECTIVE ANALYSIS OF DIGITAL TRACES 


A. REVIEW OF PRESENT IDEAS 
ie Introduction 

The mechanical bathythermograph, BT, first intro- 
duced by Rossby and Montgmory [1935] and developed by Spilhaus 
f%se] provided one of the first suecessful attempts to 
obtain a continuous graphical record of an oceanic variable, 
namely the distribution of temperature with depth. The 
expendable bathythermograph, XBT, introduced by Francis and 
Campbell [1965] increased the continuous thermal record to 
a maximum depth of 457 meters. More sophisticated devices 
with greater depth ranges have been used to record continu- 
Gusly the vertical distribution cof salinity, temperature and 
sound velocity as discussed by Magruder [1970]. Recently, 
Tucker and Shepard [1971] npadented results from continuous 
registration of light transmission versus depth. 

Continuous horizontal distribution of current speed, 
temperature, bathymetry and gravity are among the other 
variables that have been measured with sensors which are 
carried or towed in one direction at the sea surface or 
along other horizontal reference planes below or above the 


sea surface. 


Tt is important to note that although these dis- 
tributions are usually labelled as temperature (or other 
variable) versus depth (or other horizontal or vertical 


axis), the independent variable is the depth or the 
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successive positions where the sensor measures the oceano- 
graphic variable. A significant characteristic of a continu- 
ous distribution of a given variable along a length axis 
compared to time distributions is the absence of regularity 
or law of repetition, which means that distinct phenomena or 
@gifferent intensities of the same phenomenon are taking place 
along the sensed axis and makes impossible the application of 
time series analysis techniques. 

The small cost of some equipment, namely BTs, their 
reliability and easy maintenance, together with a recording 
operation that does not require great ability from an 
observer led to widespread use; an immediate consequence was 
enormous accumulation of BT traces. 

Recent development in automatic processing techniques 
of bathythermographic data is due to LaFond's [1951] 
pioneering efforts to establish a standardized method for 
processing BT traces and in predicting the need for auto- 
matic technology. Sauer [1964] described the Bathythermo- 
graph Card Processor which automatically transposes the BT 
record to an aperture card. Sauer and Hope [1967] gave a 
detailed account of the "BT Digitizer" operation which pro- 
duces semiautomatic digitization of BT aperture cards. Dale 
and Stevens [1970] presented the Fleet Numerical Weather 
Central technique to digitize the XBT traces by the 
definition of the linear portions of the curve. 

The analysis of a continuous record, BT or XBT, 


requires an acute visual sense, experience anc superior 
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understanding of oceanographic processes. The researcher 
looks for characteristic geometric features of the trace; 
usually he follows the evolution of one feature throughout 
a set of records to assure consistency of reasoning. The 
human system is slow, inaccurate for subtle variations, 
reactive to routine procedures and subjective. 

Efforts towards an automatic analysis of BT traces 
have been described among others by Boston [1966] who 
identified the most significant features in the thermal 
structure of oceanic upper layers and derived mathematical 
schemes to define the depth of the mixed layer, single and 
multiple thermoclines, transients and inversions. The 
mathematical schemes are based on: (1) an assumption of 
Gaussian or Normal cistribution of temperature as a function 
Semcepth (Figure 2), and (2) a finite difference method. 
Grosfils [1968] translated into a computer program the 
schemes above outlined and Denham [1969] applied the same 
principles to a digital analysis of internal waves. Boston 
initially developed and applied his method while examining 
nearshore internal waves. 

2.) Evaluation of Boston's Concepts 

The first part of this thesis generaiizes some of 
Boston's [1966] concepts and their presentation. An evalu- 
ation will outline and hlep to understand the purpose of the 
present research. 

a. Significant Features in the Thermal Structure 

Boston defined the significant fe>tures of the 
thermal structure obtainable with a BT within the upper 300 


meters ox the ocean as 
18 
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The Gaussian Thermocline and 
its Frequency Function [Boston, 
1966] 
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(1) Primary features or those characteristics 
common to most oceanic regions. The thermal structure is 
simply approximated by an upper isothermal layer, a thermo- 
Cline and a lower region of nearly constant temperature. 
Three parameters describe this type of structure, namely (a) 


the depth of the mixed layer (upper isothermal), Z (b) 


pod Bi 
the depth of the "bottom" of the thermocline, Zum and (c) 
the thermal magnitude of the thermocline is 


At a ae (1) 


pbs g Them as the Lemperatuxre at Zim and Tor is "the temperature 


at 4pr° 

(2) Secondary features are variations occurring 
in the primary structure at any point and classified as (a) 
EPansients or thermal fluctuations less than one degree 
Celsius and localized at the upper isothermal layer, (b) 
multiple thermoclines which are the result of heating and 
mixing conditions and (c) inversions or a type of structure 
occurring below the main thermocline where the temperature 
decreases with depth to a minimum, then increases and some- 
times regains the decreasing trend. 

b. Analytical Techniques 

The analytical methods developed by Boston are 

based on (1) a Gaussian or Normal distribution of temperature 
with depth, and (2) a finite difference approach. 


(1) The Gaussian distribution is described by 


the well-known frequency function, Figure 2, 
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where ¢ and m are constants. The function (2) has (a) 
symmetry about Z =m, (b) points of inflection at Z=mtqo 
and (c) maximum rate of change at Z =m + 5/3. Application 
of the Gaussian distribution requires that the central 
moments of odd order be zero. 

(2) The finite difference scheme implies the 
substitution of analog trace Z(T.) = f(T) by a discrete set 
of successive positions (Te Z(T_)) obtained through a con- 


Stant increment, AT, of the independent variable, Th The 


first finite difference of 2(T) will be defined as 
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and the second finite difference as 


~AZ(T) (4) 


A~z eS AZ a) 


The value of the increment AT will vary according to the 
g@e@sired accuracy of the scheme. 


e./" Application of the’ Analytical’ Techniques "ce 
the Thermal Structure 


The above concepts may be applied to the analysis 
of the thermal structure. It is assumed that the BT trace is 
replaced by a discrete set of ordered points (T,Z) from the 
surface and taken at equal thermal increment, AT. 

(1) The Gaussian Thermocline. In accordance 


with the Normal or Gaussian distribution of temperature with 
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depth, T (==), (a) the center of the thermocline is located 


VA = -m (5) 


(b) the top and the bottom of the thermocline are considered 
to be at the points of the maximum rate of change of the 


slope of the T (2) curve or at 


pe ee (6) 


N 
I 


Za 3m =G Yo (7) 


(2) The Non-Gaussian Thermocline. The analysis 
of the non-Gaussian thermocline makes use of the variation 
in sign of the second finite difference scheme, equation (4), 
EeeeinGg the point off inflection of the (T,2)\curve. This 
Perme tdentifies the center of the thermocline, 2 .. The cop 


of the thermocline, Z is determined using the technique 


oul ied 
outlined in (1) and generating an artificial Gaussian dis- 
tribution around Z op by reflection of the upper part of eEhe 
trace (Z =2ZT). The bottom of the thermocline, Zh n is 
obtained in the same way but making use of the lower part 
Of the (1T,Z) curve, (2 = Z_.). Figure 3 shows the successive 
steps of this analysis. 

(3) Transients are recognized by the sign changes 
G@retne second finite Gueeerence m7 (1): A sign change from 
positive to negative signifies the end of one transient and 
ene beginning of the next one: 


(4) Multiple Thermoclines. An examination of 


the sign changes of the second finite differe:ce, as made in 
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Figure 3. A Non-Gaussian Thermocline 
Showing Top, Center, and 
Bottom as Determined by 
Finite Difference Procedure 
[Boston, 1966] 
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(2), will reveal the points of inflection of the (1T,Z) curve. 
These correspond to the centers of the thermoclines. The 
top and the bottom of each thermocline are found applying 
the method explained in (2), The separation between thermo- 
Cclines is determined by the location of successive points 
where the first finite difference attains a maximum value. 
Examples of multiple thermoclines are illustrated in Figure 
4. 

(5) Inversions. The sign change of the first 
finite difference of the (T,Z) curve indicates the existence 
of an inversion. The upper and lower limits of any inversion 
Aadineste a continuum where the representative trace is 
double-valued (same temperature at two depths) but outside 
the region the distribution of temperature with depth is 
single-valued (one temperature and one depth). The main 
features investigated by this method are shown in Figure 5. 

d. Remarks about Boston's Concepts 

(1) The problem of automatic processing of 
bathythermograph data is well stated and represents a 
remarkable contribution. 

(2) The objective digital analysis of the thermal 
structure requires application of several methods, each one 
is independent from the others and each describes a distinct 
thermal feature. More general objective schemes to analyze 
other features in the thermal structure would be even more 
useful and might be extended to the objective study of other 


oceanographic variables with vertical or horizontal axial 


aGistributions. 
24 
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TABLE IV. 


Multiple thermocline analysis, 
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Figure 4. Thermal Structure with One, Two, and 
Three Thermoclines [Boston, 1966] 


‘ 
; : 
7 i ob beg 
= 4 ‘ : Ai 
oF 4 
' : a 3p eRe Th: 
4 Vern rrr ry 
i A a - 
1 4 
: an on y id 
? oi by : \ 
' . 
{ ) y 
- : 
Lat LMT Bat ieee ‘ 
= ' 
a a9 Mir eegou 7 
; moms st _ 4 w 
~ Sie ee 
7 \ 
a | 5 
>I ba 
— : — \ 
f ea aT 
A : \ 


/ . as E 
hee. r . Rocha thes hachnlnahoel Qi 


= : _ a eee? Jets 


i's 
f 


DEPTH (METERS) 


40 


80 


120 


160 


200 


240 


260 


Tmin 
Zc 
UPPER BOUNDARY 
Z(Tmin) FOF 
DEPTH OF (MAGNITUDE ) 
Tmin 
LOWER BOUNDARY 
Zc2 
5 10 1 20 


TEMPERATURE (°C) 


FIGURE 5. Inversion Illustrating Main 
Features [Boston, 1966] 
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(3) The Gaussian or normal distribution Be 
temperature with depth constitutes an interesting investi- 
gation... Its ,application to the definition of .a .thexmo- 
cline is restricted depending upon whether or not a 
condition of symmetry around the first moment does exist. 
If this provision is not fulfilled there is not any thermal 
normal distribution with depth. However, in this case of 
non-Gaussian thermocline the top and the bottom of the 
thermocline are determined by means of an artificial 
reflected construction of normal distribution around the 
center which distorts slightly the initial information. 

(4) The finite difference technique explores 
conveniently the behavior of first and second order schemes 
which give a familiar feeling for fluctuations of slope and 
Bones of inflection of the original trace., The development 
and study of higher order schemes may indicate interesting 
variations associated with significant features of any 
curve. 

(5) The digitization of the analog trace by the 
selection of depths at uniform temperature intervals 
appears to be advantageous because the data is ina form 
where it may readily be processed by the finite difference 
technigue. The thermal increment should be chosen as small 


as possible but compatible with the accuracy of, the Sensor. 


B. AN OBJECTIVE MODEL 
The primary purpose of this chapter is the presentation 


and development of general concepts which ma} be utilized 
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in the objective analysis of any given two-dimensional 
trace. 
Lt. Significant Features in a Two-dimensional’ Trace 
A two-dimensional trace presents the distribution of 
any oceanographic variable, V, with depth or along a hori- 
Zontal axis, S, aS shown in Figure 6. A pecularity of the 
Erace Consists in the continuous increase in’S with 
Eluctuations in V. Such curves (V,S) show three general 
characteristics, namely: (1) trend, (2) varicline, and (3) 
inversion. 
a.! Trend 
The trend indicates the general direction of 
development of the curve (V,S) in Figure 6: (1) between the 
beoznning of the record, a, and q, the variable Vv is 
decreasing at different rates to the left of the figure; 
(2) 9 -— h; the variable V is increaSing and the trace is 
directed to the right; and (3) the final position of the 
curve between h and i again shows a decrease of V. 
bi Varicline 
Varicline was derived from the comoosition of 
the words: relative "variation" of the slope or "inclination"; 
it identifies portions of the curve where the slope increases 
relative to the neighbor regions. In Figure 6, parts b- Cc 
and d - e of the curve are variclines. Considering vari- 


cline ad - e, the limiting parameter values are (ir) the top, 


b, 


= (8) 
Siy S(d) 
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Figure 6. A Two-dimensional Trace 
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(2) the bottom, e, 


Shy =a Set) (o) 

(3) the range, AS, 
ASy = S(e) - S{d) (10) 

(4) the center, 

Soy  $(a) tAS/2 (1) 
(5) the variable magnitude, AV, 

AV, = vie) - Vv(d) (12) 
and (6) the gradient, 

Gy = AV, /AS,, (13) 


ec. ~ inversion 

An inversion identifies a segment of the curve 
(V,S} where there exists an approvimately symmetrical dis- 
tribution of the oceanographic variable, V with length S 
around a point, namely location g in Figure 6. The point of 
aymmeury 15 found by examination of the variation’ of the 
trend: it is a relative maximum or minimum of V. Referring 
EOeragure 6 again, the point g constitutes a relative mini- 
mum. The most important parameters of an inversion are (ar) 


Ene cop, £, 
S(£) (14) 


op) 
MN 


(2) the bottom, h, 
S(h) (15) 


Sr 


(3) the location of the extreme value of V in the inversion, 


g, 


il 


16) 
Sor S(g) ( 
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(4) the variable magnitude, AV 


iE 
Av, =|v(g) - v(£)| (17) 
=lv(g) - v(h)| (18) 
and (5) the range, AS, 
AS. - s(ilie ste! (19) 


Ze ane Pinite:r Difference Technique 


The finite’ difference technique investigates the 
Variation of a two-dimensional distribution. The original 
analog record is previously digitized according to semi- 
Seutomatic *procedures. The digital trace consitutes an 
ordered and discrete set of positions. 

a. The Preparation of the Mesh 

The finite disference technique requixes) an 
organized two-dimensional distribution or mesh such that 
the values of one variable V, chosen as independent, are 
read at constant intervalav together with the correspondent 
values of the other Saisie S(V), dependent, Figure 7. The 
magnitude of the incrementAV is directly related to the 
accuracy of the record and it is unrealistic to select a 
small but incompatible value. The digital trace may present 
the structure as a satisfactory mesh with values at constant 
intervals; if not, it must be altered accordingly. Assuming 
that the digital trace has a linear variation, the mesh 
points will be obtained by a simple linear interpolation. 
However, the selected constant interval may introduce an 


alteration to the original trace and AV must be small enough 


as shown in Figure 8. 
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The Linear Interpolation of a 
Digital Trace: “A.” Digital DTrace, 
B. Mesh (increment V = 2 arbitrary 
untes)2'"C2" Mesh (inerement V7" 


arbitrary unit) 


Figure 8. 
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bd. Finite Ditference Schemes 

A mesh consists of an ordered set of points at 
fixed intervals of variable and two consecutive positions define 
one linear segment. Each segment is specified by (a) the 
length component and (b) the variable component as illustrated 
an Figure 9A. 

(1) First Finite Difference of V 

The first finite Gireerenee of V(V_) abe 

defined as 


EMA) = VV = Wile) (20) 


n+1)? 
and may be positive (pointing to the right) or negative 
(directed to the left), as shown in Figure 9C. 
(2) First-Finite Difference of S 
The tirst finite difference of S(V_) is 
stated as 


AS(V_) = SV = S(V_) (21) 


n+1)? 
and is always positive, Figure SBE 
(3) Second Finite Difference of S 
The second finite difference of S(v_) is 
described as 
A*s(v,) =As(v 


ae (22a) 


= Sy 25(V4,) * Sv) (22b) 


nt2) 
and employs three points, n, ntl, nt2, from the mesh. The 
three points define two linear segments, whose relative 
arrangement will be related to the negative, positive or 


null value of a*s(v.), Figure 10. 
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Figure 9: 


The Mesh and Its Components: 
A. Segments of the Mesh 

B. Length Components 

C. Variable Components 
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Figure 10. Schemes of the Second Finite 
DrLLerences: A. ‘Constant 
Slope, B. “Increasing Slope, 
C. Decreasing Slope 
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(4) Third Finite Difference of S 
The third finite difference of S(V_) is 
specified as 


A*s(v.) =A7s(v_4,eoaesiyD (23a) 
= ASV 25) = o2he(i a 1) ee eeey 


= 5S (a, Je Sony - 3S(V 


n+2) ae 


‘i Sie) 


(23¢) 


and requires four points, n, ntl, n*2,-nt3, £rom the mesh. 
The As(V.) summarizes the relative disposition between 
EnLee PTinear Segments, defined by fourspoeints.* Five basic 
schemes of the A*s(v_) are presented in Figure ll. 
(5) Higher Order Finite Difference of S 
Higher order finite difference of S, 


malt eon ee ) (24) 


k _ » 
- S(v.) = nel n 


S(V 
will demand more points from the mesh, namely the number of 
the order plus one and concurrently the analysis of the 
relative structure will increase in complexity. The develop- 
ment of the finite difference schemes up to the third order 
seems to constitute adequate advancement for the operational 
meeas in the present research. 
3. Two-dimensional Objective Analysis 

The finite difference schemes above introduced, will 
be applied to the significant features of a two-dimensional 
trace and illustrated by schematic diagrams. | 

a. Lrend 


As explained before, the trend shows the tendency 


of the independent variable, V, which is investigated by the 
Si 
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Schemes of the Third Finite Differences: 

A. Constant Slope, B. . Successive Increase 
of the Slope, C. Successive Decrease oF 
the Slope, D. Successively Decreasing and 
Increasing the Slope, E. Successively 
Increasing and Decreasing the Slope 
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use of its first finite difference, AV(V_) expressed in (20). 
Assuming that V increases to the right, as shown in Figure 
6, the analysis of V may be summarized as follows: 
(1) Variation 
(a) AV(V_) On: Vera < V., V is decreasing 


n 
or directed to the left of the graph. 


(b) Av(V_) = Oe Vaiyes Vin! V is constant 
or parallel to the length axis. 

(c) AV(V,) Sr One Viet > V,, V is increasing 
eretending to -the-right. 

(2) Extreme values 

The change inesigniof thecfirst fLimite 
difference of V indicates that the trend has changed direction 
end defines an extreme feature: 


(a) ecViiv < O and Vin) > 0 : relative 


n+1) 
maximum at point ntl. 

(b) Vi eh) > O and WV a) < 0 : relative 
minimum at point ntl. 

(c) The largest and the least of the V values 
above found with the values of V at the beginning and end of 
the trace give the absolute maximum and minimum of the trace. 

b. a Vabicline 
The basic concept of the objective varicline is 
presented in Figure 11E composed of three segments and, 
although their relative arrangement is unchangeable the 


middle part can be divided into more than one segment. The 


analysis of the characteristics, 
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(1) a*v(v,) < 0 : top, of varieline at, point. ntl. 


(2) Av (Vv ) > O : base or bottom of varicline 


ntl] 
Ae point nt2 is indicated in Figure. 11; 
Ge Tnversion 

The analysis of the trend indicates whether or 
not there are maxima and minima present in the record. The 
maxima and minima (excluding the values of V in the beginning 
and the end of the trace) evidence areas of approximate 
symmetry of V around the extrema. Referring to Figure 6 
the analysis of the trend gave (a) the absolute minimum at 
g, (b) a relative maximum at h, and (c) the absolute maximum 
value in V was found at a, the beginning of the trace. It 
is evident that locations £ and h define a region of approxi- 
mate symmetry of V around the relative minimum g. An inver- 
Sion is characterized by a portion of the (V,S) curve where 
Ene trenG changes sign, illustrated by the layer g = h-in 
Figure 6. The limits of the approximate symmetry are 
established by comparison of the variable magnitudes, AV'S 
(17 or 18) determined using the point cf symmetry of the 
extrema in consideration, g, and the adjacent positions 
defining extreme values in V, namely a and h. The least AV 
is chosen as the variable magnitude of the inversion. 

ade. "eonplex type 

The complex type of trace combines the basic 
investigations of trend, variclines and inversions. The 
method of procedure is outlined in three stages: 


(1) Objective analysis of the trend. 
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Figure 12. 


TEMPERATURE —— 
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4 


Computer Output, Plot of a Digital Trace 
Showing the Tops (Squares) and the Bottoms 
(triangles) of Analyzed Variclines. (The 
Control of the conditions limiting the 
bottom of a varicline allows: A. jdentz— 
fication of all’ the variclines present, and 
B. definition of the significant features) 
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(2) Objective analysis of inversions. 
(3) Objective analysis of variclines, which may 
be conducted throughout the entire trace or be restricted 
to those portions outside the symmetrical regions. 
4. Tri-dimensional Objective Analysis 

The two-dimensional objective analysis studies the 
distribution of a given variable V along a length axis S. 
Mne availability of a collection of traces spread out in 
eeme, DD, or along,a horizontal, L, ox vertical, S, direction 
illustrated in Figurel3, makes possible the investigation of 
a tri-dimensional distribution of V in (a) length (V,S,L) 
and in (b) length and time (V,S,D). An objective analysis 
of this tri-dimensional distribution is possible if the law 
Sf Variation of V is Known cr the interval fin time--.o3 
length) between consecutive traces is small enough to assume 
Seeanear interpolation or digitization of the new trace: 
One of the above requirements is fundamental to generate a 
mesh and to apply the concepts of the two-dimensional 
objective analysis to each one of the new traces, shown in 
Figure 13B. 

ae The weoblemioft Filtering 

Usually a trace shows a complicated form with small 
Variations superimposed on a general pattern. The small 
fluctuations may have physical significance or may simply 
represent an artificial response of the sensor; however, 
it may be impossible to distinguish whether a regular 


feature is valid or invalid. The most common 2pproach 
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Tri-dimensional Analysis: A. The 
Tri-dimensional Distribution, B. 

A New Trace, C. The Projection of 
Two-dimensional Traces Along a Given 


Axis 
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Figure 14. 
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The Filtering Parameters of Variclines 
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Fagure-152°¢° Computer Output Plot ‘of “the Linear 
Field of Isotherms for One Cruise or 
the Average Linear Field of Isotherms 
for a Given Number of Cruises 
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consists in the identification of the general pattern of 
the trace and this outlines a problem of scale: if the 
investigator is interested in microstructure, his main 
efforts are concentrated in the improvement of the accuracy 
Of the sensor. The technique of suppressing the small 
features has not found definite and complete solution. It 
reveals a compromise established between the need to avoid 
small fluctuations and the unfortunate results when those 
variations were physically meaningful and their filtering 
has introduced artificial values. The best approach requires 
a good understanding of the geophysical phenomena and a 
Suitable application of the filtering methods. 

in this thesis study, analytical filters (such as (the 
binomial expansion and the least Square polynomial apprexi- 
mation were used to smooth the generated mesh but they 
introduced severe modifications which (as presented in 
Figure 16B) completely deformed the original information, 
and their further application was avoided. 

The discussion of the filtering capability of the 
finite difference technique constitutes the purpose of the 
Present chapter. 

a. The Mesh 

The generated mesh depends upon the selection of 
the constant incrementaVv: (a) a large value erases the small 
features but introduces distortion in the digital trace, and 
(b) a small value (compatible with the accuracy of the sensor) 
keeps the details and reproduces closely the di vital trace, 


as shown in Figure 8. 
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Figure 16. Computer Output Plots of: A. The 
Interpolated Digital Traces Along a 
Cross Section fora Given CEulse, 
and B. ‘Ihe Fifth Order Lease 
Square Polynomial Approximation of 
the Same Digital Traces 
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b. Significant, Features 

Generically, Significant, features like inversions 
and variclines may be filtered out by specification of 
filtering parameters or the minimum values of (a) the 
defining parameters (8 - 11) and (b) the separation in 
length between consecutive and identical features (bottom of 
one and the top of the next). 

Specifically, the filtering of variclines may be 
obtained by: 

(1) The conditions limiting the top and the 
bottom of variclines (especially the bottom) and expressed 
ett. b.35.D., May be modified and extended to more points 
of the mesh above and below the feature resulting in its 
more clear definition, as evidenced in Figure 12. 

(2) The filtering parameters illustrated in 
Figure 14: 

(a) ALPHA or the minimum variable magnitude 
mimtcing the varicline. 

(b) BETA or the minimum gradient defining 
eae Varicline. 

(c) DELTA or the minimum separation in 
length between consecutive variclines. 

(d) THETA or the minimum gradient selecting 
the variclines to be merged, using simultaneously the SIGMA 
parameter. 

(e) SIGMA or the maximum separation in length 
between variclines to be merged and already defined by the 


THETA parameter. 
48 
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C. A COMPUTER PROGRAM 

The investigation of the verticaly;diseEri bution) ef tempera— 
ture in the upper oceanic layers can provide a digital model 
which applies the concepts of objective analysis described 
above. 

Usually the data available is too non-uniformly distributed 
in time and length to make possible the use of a tri-dimensional 
objective analysis; however, a Eyosaieneserer type may be 
weplized to examine each trace: and to: plotsthe resultwetethe 
investigation (of all the traces obtained in a cruise or 
passage along a cross section) in a way which allows a 
subjective analysis. 

Anointroduction Of the model, its characteristics: and 
Seeration is summarized in the next presentation. 

fonehbescription 

a. The Basic Computer Model 

The computer model is organized in two parts: 
(a) the main program, and (b) the subroutines. The main 
program presents in a logical way a sequence of events con- 
cerning the execution of the job. This arrangement is a 
compromise between the space available in the computer 
memory and the least time of execution. The subroutines 
deal with specific stages of that execution. 

The following description outlines the major 
operations of the model and their implementation must be 
monitored closely by the use of the flow diagram presented 


in Appendix B, with the symbols explained in / .:pendix A. 
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(1) “START toc" NORDER 2=au! 

The computer program receives the general 
information about the number identifying the cross section 
in study (ICSEC), the number of cruises per cross section 
(ICRUI), the number (ISO) and the values (TI) of isotherms, 
the mesh interval (DTW) and the number (NBT) and location 
(X) of stations for the cross section. 

(2) "DO 460" Loop 

The computer program processes the complete 
set of information about the thermal structure and wind field 
Bom a Given Cress Section; cruise after cruise: 

| (a) "ITITLE'S" reads the titles of the plots 
drawn per cruise. 

(b) "DO 255" loop reads the location (ISTA), 
the number (NOP) and pairs of points (ZO,TO) defining each 
Gigital trace in the cruise and processes the’ information 
about the depth of each isothermal value (GTVALU or ARTFT1). 
This scheme applies to all digital traces of a cruise. 

(c) "CALL PLOTIT" receives the above gener- 
ated data for the whole cruise and plots the field of 
isotherms (TI). 

(d) "DO 350" loop reads the location (ISTA), 
the number (NOP) and pairs of points (Z0,TO) defining each 
digital trace in the cruise and prepares the mesh for that 
trace (LINTPS), plots the obtained mesh (PLOTBT) and analyzes 
it to identify thermoclines (THERMC). 

(e) "CALL PLOTFT" receives the above generated 


data for the whole cruise and plots the field of thermocline 


isopleths. 
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(f) "DO 450" loop reads the location (ISTA) 
and the wind observations (WIND); determines the components 
of the wind along and perpendicular to the next coastline 
at the day of XBT/BT drop and one day before (WINDO1), two 
days before (WINDO2), three days before (WINDO3), five days 
before (WINDO5) and ten days before (WIND10); plots the wind 
and its components. 

(g) "NPLOT = 3" to "460" completes the draw- 
ing of the wind information for a cruise and sets in the 
process for a new cruise. 

(3) PEKONTR! Oto. TKONTRES OY 
This step controls the execution of the 


Pomeshing partoof the program: «Lf the ZKONGR.value is zero 


(4) SEREODr oink." otol "CALL. STATESS 

The computer program reads the required 
information to start the statistical study (STATCS) between 
any two parameters (KEK = 1, KEK = 2) and to plot the scatter 
diagram and the least square regression lines in the same 
Srapht {ITITLS).. ~The continuation of the statistical study 
relies on the value of KRK read if different from zero. 

(6). “LONTR" to "TLONER = 0” 

This part controls the execution of the-last 
part of the program. If the ILONTR value is zero, the job 
ends. 

(6) VE TTTh6” oot PLIMaE" 

The program receives the gene: ated data of the 

whole cross section and plots the average field of isotherms 


(PLTMIT) with a printed label (ITITL6). 
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(7%); XECSEC! sto. SSTORAGE” 

The program indicates that the study of jthe 
cross section (ICSEC) is completed and ends the job by the 
storage in a disk of the original and generated information 
(STORAGE). 

b. Other Subroutines Available 
Other subroutines are available to improve or to 
modify the operation of the basic modetaiae keeping its main 
SEructure.. 
(1) Subroutine TREND 

The subroutine TREND, Appendixes D,-G. and M, 
analyzes the digital trace using the finite difference 
technique to locate the maxima and the minima. 

(2) Subroutine INVERS 

The subroutine INVERS, Appendixes E, G and M, 
determines the inversion parameters for each minimum 
identified by subroutine TREND. 

(3) Subroutine PROPAR 

The subroutine PROPAR, Appendixes H and K, 
computes the friction velocity, roughness parameter and 
Stability length from wind speed and air-sea temperature 
difference data. 

(4) Subroutine ITRATE 

The subroutine ITRATE, Appendixes I and L, 
calculates the friction velocity, inflow angle, roughness 
parameter and stability length from the large scale synoptic 


parameters. 
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(5) Subroutine DRAW 
The subroutine DRAW, Appendix J, neutralizes 
the "IBM subroutine DRAW" avoiding output plots. 
Z2.0* Input 
The form of input to this program consists in data 
deck. 
a. Main Program 
The input to the main program is identified by 
the symbol READ shown in Appendix A. The flow diagram of 
the main program, Appendix B, gives the ordered sequence of 
the> data deck. The format of individual input is specified 
in the computer program, Appendix F. 
b. Subroutine THERMC 
The subroutine THERMC requires the insertion of 
chosen values for the filtering parameters (ALPHA, BETA, 
DELTA, SIGMA, THETA) before running the program. 
c. Subroutine PROPAR 
The subroutine PROPAR requires special data 
described in Appendix H, which must be read together with the 
wind information as shown in Appendix kK. 
d. "Subroutine ITRATE 
The subroutine ITRATE needs large scale synoptic 
information pictured in Appendix I, which must be read 
together with the wind data as presented in Appendix L. 
a. ~OuEDUE 
The forms of putput from this program are (1) 


Beintouts, (2) plots, and (3) storage in disk. 
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a. Printouts 
The printouts are identified by the symbol 
PRINTOUT, in Appendix A. The flow diagrams show the occur- 
rence of printouts and their formats are specified in the 
computer programs. 
The inclusion of the DUMMY card, Appendix P.a. 
Meutralizes this form of output. 
b.” PLoOts 
The basic program presents six types of plots 
which appear according to the execution of the following 
subroutines: 
(1) 2 PLOLET 
The subroutine PLOTIT plots the finear freig 
Of isotherms per cruise, as shown in Figure 15. 
(2) PLOTBT 
The subroutine PLOTBT plots the interpolated 
Gigital traces obtained in a given cruise, Figure 16. 
(3) PLOTFT 
The subroutine PLOTFT plots the linear field 
of isopleths defining the top, the mean depth and the bottom 
of the thermoclines analyzed from the interpolated digital 
Eraces Of a given cruise, Figure 17- 
(4) PLTMIT 
The subroutine PLTMIT plots the average 
linear field of isotherms for a specified number of cruises, 


Figure 15. 
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Figure 17. 


CROSS SECTION AXIS —> 


Computer Output Plots of the Linear Field 

of Isopleths Defining the Top (T), the 

Mean Depth (M), and the Bottom (B) of 

the Three Shallowest Thermoclines for 
Different Values of the Filtering Parameters 
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(5) RLINE and STATCS 
The’ subroutines RLINE and STATCS plot the 
least square regression lines (RLINE) and the scatter dia- 
gram (STATCS) in the same graph to illustrate the statistical 
relation between two given variables, Figure 18. 
(6) WINDOO, WINDO1, WINDO2, WINDO3, WINDO5, WIND10 
The subroutines WINDOO, WINDO1, WINDO2, 
WINDO3, WINDO5 and WIND1O plot the wind vector and its com- 
ponents (along and perpendicular to the direction of the 
near coastline) measured at given station and day of XBT 
observation and at one, two, three, five and ten days before. 
The drawings are made in the same diagram for a specific 
Puce, Figures 19. 
(7) DRAW 
The insertion of the subroutine DRAW, Appendix 
weowill avoid this type of output. 
Ge Storage in Disk 
The initial and processed data are stored in 
allocated space of disk depending upon the inclusion of 
the Job Control Language instructions shown in Appendix O.b, 
e,a. 
4. Operation 


a. Characteristics of the Basie Program 
(1) Source Language 
The computer program is written in FORTRAN 
Iv (H) and the calculations are performed in real single 


precision arithmetic. 
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Figure 19. Computer Output Plot of the Wind vector 
and its Components (along and perpendicular 
to the direction of the near coastline) 
Measured at Same Station and Day of XBT 
Observation and One, Two, Three, Five 
and Ten Days Before 
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(2) Subroutines and Functions Subprograms 


The computer program uses subroutines and 
functions subprograms that had been compiled or been added 
to a library under OS/MVT. If the source deck of subroutine 
MISR had not been precompiled, it might be accessed (in 
Naval Postgraduate School Computer Facility) by concatenation 
in the FORTRAN compile step using the Job Control Language 
cards shown in Appendix N.d,e. | 

(3) Core 

The basic computer program processes infor- 
mation for twenty cruises, each containing twelve stations 
and requires 288 K in core dimension, demanding the FORTRAN 
compilation step 150 K which must be requested by means of 
Ene Job Control Language card presented in Appendix N.b. 

The space available is divided in five 
"COMMON" labelled areas which enables (a) the exchange of 
information between the main program, and (b) an economic 
and profitable use of the computer memory. 

(4) Time 

The time required to process (without 
Statistical study) a set of seventeen cruises, each with 
nine stations, is about three minutes and thirty-four 
seconds. This time is distributed by (a) the compilation 
(one minute and forty-one seconds), (b) the linkage (three 
seconds), and (c) the execution (one minute and fifty 


seconds). 
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(5) Disk 
The storage of the initial and processed 
fnformation in a disk is obtained by the use o£ Ene Dob 
Control Language cards shown in Appendix O.b,c,d. The 
requested space in disk corresponds to one cylinder or IBM 
2314. 
b. Equipment 
The basic equipment used was: 
(1) IBM 360/Model 67, OS - MVT with a 


(2) CALCOMP Plotter, Model 765. 
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III. APPLICATION OF THE OBJECTIVE MODEL 


The computer model was tested by the use of bathythermo- 
graphic and wind information available from Fleet Numerical 
Weather Central for the geographical area shown in Figure l. 

The XBTs were taken during three years of passages 
between Canary Islands and Ferrol (Bay of Biscay) and 
Canary Islands and Barcelona (Mediterranean Sea) by Spanish 
Ships. The XBT traces, in digital form and FORTRAN. format 
compatible with this program, were organized in three cross 
Bections, Gach one describing the normal routing of the 
Survey ships and containing ten stations. In total, about 
750 XBT traces were ready for further processing. Values 
for the filtering parameters of the thermocline analysis 
were previously chosen, namely 

ALPHA = 0.3°C 

BETA = 0.015°C/meter 

DELTA —- 10.0 meters 

THETA = 0.027°C/meter 

SIGMA = 25.0 meters 
and applicable to all seasons and stations of the area. 
Finally, the objective analysis of the bathythermograph 
information was performed. 

For each station of the cross sections, geostrophic wind 
data were computed from surface pressure analysis in FNWC. 
The wind obtained was not corrected for instability due to 


lack of information concerning the air-sea temperature 
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difference and according to Cardone [1969], such correction 
is fundamental to assure reliable wind data, specifically 
periction velocity”. Without this transformation, further 
handling of wind was decided to be inappropriate. 

The final part of the application of the model, which is 
a statistical study between the wind and the thermal fields, 
as well as the computation of the first, second, third and 
fourth statistical moments (to infer isnot statistical 
distributions) of each field for a given station of the 
area could not be physically concluded due to the lack of 
proper wind data. Otherwise, the difficulties pointed out 
the importance of valid wind information and motivated the 
inclusion of subroutines, developed by Cardone [1969], to 
Saiculate "friction velocity" from large scale synoptic 


information (ITRATE) or ship reports (PROPAR). 
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IV. CONCLUSIONS 


A FORTRAN IV program was developed and tested to handle 
large amounts of wind and bathythermographic information in 
an automatic way without distorting the initial data and 
restricting the human participation to a preparation phase 
dealing with organization of the information available and 
selection of values for the parameters of the model. The 
program performs objective analysis of XBT/BT digital traces, 
processes wind information, studies the statistical character- 
fetees Of the original and processed data and plots the 
generated fields. 

The objective analysis of XBT/BT traces makes use of 
finite difference schemes which requires preceeding trans- 
formation of the analog trace into digital form. The same 
concepts may be extended to analyze the distribution of any 
environmental variable (similar to the bathythermograph) 


along vertical or horizontal directions. 
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V. RECOMMENDATIONS 


It is recommended that bathythermograph data be analyzed 
by an objective method before any subjective interpretations 
are made. 

The objective analysis of XBT/BT traces may be performed 
m@board ships by the use of digital computers requiring a 
small memory but, the pre-digitization of the analog trace 


constitutes a fundamental prerequisite for processing. 
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APPENDIX A 


FLOW DIAGRAM SYMBOLS 


JoUcoVOouud 


LIMyTS 


READ 
PRINT OUT 
DECISION 


ately 


STORAGE.IN DBiSts 


CONNECTION 
SUBROUTINE OR FUNCTION 
SIN a} jolly je jalale 


FPRECESS 
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APPENDIX B 


FLOW DIAGRAM OF THE FORTRAN 
IV BASIC COMPUTER PROGRAM 
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APPENDIX C 


FLOW DIAGRAM OF THE THERMOCLINE ANALYSIS 
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APPENDIX K 


MODIFICATION TO THE FORTRAN IV BASIC 
COMPUTER PROGRAM (APPENDIX F) TO INCLUDE 
THE SUBPROGRAM PROPAR (APPENDIX I) 
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MODIFICATION TO THE FORTRAN IV BASIC 
COMPUTER PROGRAM (APPENDIX F) TO INCLUDE 
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JOB CONTROL LANGUAGE (JCL) FOR 
THE HEAD PAR’ OF THE FORTRAN IV 
BASIC COMPUTER PROGRAM (APPENDIX F) 
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